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Abstract 

There is strong evidence in favor for zonal flow suppression of the 
Ion- Temperature-Gradient (ITG) mode turbulence, specifically close 
to the linear stability threshold. The present letter attempts to an- 
alytically calculate the effects of zonal flow suppression of the ITG 
turbulence through deriving a modified dispersion relation including 
the back-reaction of the zonal flows on the ITG turbulence based on 
the quasilinear theory. The results are manifested in a reduction of the 
linear growth rate and an increase in the effective linear ITG thresh- 
old. 
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I Introduction 



Anomalous transport remains one of the main concerns in magnetically con- 
fined plasmas. The anomalous transport in the core is commonly attributed 
to Ion- Temperature-Gradient (ITG) mode turbulence, however, other modes 
can not be totally excluded [I]. An important critical concept is the trans- 
port regulation and transport barrier formation by zonal flows. The zonal 
flows are radially localized (radial wave number q x ^ 0) and strongly sheared 
flows in the poloidal direction (poloidal and parallel wave number q y = 0, 
Q\\ = 0) [2]- [6]. The large scale zonal flows are believed to be excited by the 
small scale fluctuating drift wave turbulence via Reynolds stress. 

Several models for zonal flow generation from turbulence have been pro- 
posed, wave kinetic equation models, reductive perturbation methods and 
coherent mode coupling models [6]. The coherent mode coupling model used 
in this particular study has been adopted and fitted to this field of study 
through the recent years [6]- [13] . It has been shown that the model is in 
good agreement with the corresponding model derived using the wave kinetic 
equation [T3| and the model has successfully contributed to increased under- 
standing of the weak generation of zonal flows in recent measurements at the 
Columbia Linear Machine [14] . 

In this letter the effect of the back- reaction of zonal flows on the ITG mode 
turbulence is analytically investigated. Previous models on back-reaction of 
zonal flows on turbulence have mainly been concerned with linear estimates 
of shearing [6], [TSj. Dynamical models for the coupled system of zonal flows 
and turbulence have also been proposed [6], [TJ, however, no explicit estimate 
of the reduction of linear growth rate or up-shift of the critical gradient (e.g. 
VTj) due to zonal flow stabilization have been made [TB]. This letter is mainly 
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concerned with estimating quantitative values of reduced linear growth rate 
and of the up-shift in the critical gradient. The inspiration for the research in 
the present letter is from the study of self-organization of plasma turbulence 
where zonal flows inhibits anomalous diffusion [3] and the phenomenon of up- 
shift in the threshold value for the critical gradient or Dimits shift that has 
been observed in many non-linear simulations of ITG mode turbulence [16J. 
A possible explanation for the Dimits shift is also recognized as one of the 
main yet unsolved physics issues for the understanding of zonal flows [T7] . 

Based on the principle of quasilinear theory, the strategy of the letter is to 
first consider the zonal flow generation from the ITG background turbulence 
and then consider a small change in the drift wave perturbation calculated 
through the non-linearities. This small non-linear correction term is found 
by letting the zonal flow couple back to the side-bands resulting in a change 
in the background drift wave system. This back-reaction term can be inter- 
preted as a change in the threshold of the critical gradient (rj ithNL ) for the 
drift wave system. Moreover, the term changing the background is shown to 
cause a decreasing growth rate and an effective up-shift of the threshold. 

The letter is organized as follows. In Sec. II the model for the ITG 
mode is presented together with the model of zonal flow generation and 
back-reaction. In Sec. Ill the zonal flow saturation level is estimated and 
the results are presented in Sec IV. We provide a conclusion in Sec. VI. 

II Governing equations and model 

The ITG mode turbulence is modeled using the continuity and temperature 
equation for the ions and considering the electrons to be Boltzmann dis- 
tributed and the quasi neutrality is used to close the system [18]- [19] that 
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has been successful in reproducing both experimental [20J and non-linear 
gyro- kinetic results [I]. The advanced fluid model is expected to give quali- 
tatively and quantitatively accurate results in tokamak core plasmas or the 
flat density regime. For instance, this was shown in the plot of the slope of 
the ITG threshold as a function of e n in the flat density regime (see Figure 1 
in Ref. [21] or Figure 4 in Ref. |22j) where this model recovers the full linear 
kinetic result within 5%, far better than the prediction from an approximate 
kinetic model using a constant energy approximation in the magnetic drift. 
In the present work the effects of parallel ion motion, magnetic shear, trapped 
particles and finite beta on the ITG modes are neglected. It has been found 
earlier that the effect of parallel ion motion on the ITG mode is rather weak, 
however, it has recently been shown that magnetic shear and parallel ion 
motion may modify the zonal flow generation significantly |23j. The zonal 
flow generation from Trapped Electron Modes [21] and the effects of parallel 
ion motion [H] on zonal flow generation have been treated earlier. 
The continuity and temperature equations are, 
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The system is closed by using quasi-neutrality with Boltzmann distributed 
electrons. Here [A, B] = (dA/dx)(dB/dy) — (dA/ dy)(dB / dx) is the Pois- 
son bracket. With the additional definitions n = 5n/no, <fi = eS<p/T e , 
Ti = 5Ti/TiQ as the normalized ion particle density, the electrostatic po- 
tential and the ion temperature, respectively. In the forthcoming equations 
r = Ti/T e , = p s c s y/L nj p s = c s /Vt ci where c s = JT e /rrii, Vt ci = eB/rriiC. 



We also define Lf = — (dlnf / dr)" 1 , rji = L n /L T ., e n = 2L n / R where R is 
the major radius and a, = r (1 + r)i). The perturbed variables are normal- 
ized with the additional definitions h = [L n j p s )Sn/no, (ft = (L n /p s e)S<p/T e , 
Ti = (L n / p s )STi/T i0 as the normalized ion particle density, the electro- 
static potential and the ion temperature, respectively. The perpendicular 
length scale and time are normalized to p s and L n /c s , respectively. The 
geometrical quantities are calculated in the strong ballooning limit (9 = 0, 
Qi (9 = 0, ac) = 1/ac where & (9) is defined by u D (9) = u+e n gi (9)) 
The approximate linear solutions to Eqs 1 and 2 are, 
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Finite Larmor Radius (FLR) effects in the rjuh are neglected, here, however, 
in the remaining part of the letter the linear state is solved numerically and 
FLR effects are included. 

Following the same procedure as in Refs [Tj, [13]- [14] by choosing a 
monochromatic pump wave that generates a secondary convective cell. This 
gives a problem with four coupled waves where the zonal flow seed may be 
unstable and grow exponentially. Here, the sidebands may couple to the 
pump wave to give the zonal flow component. No Boltzmann distribution 
for electrons is given to derive the third order dispersion relation for the zonal 
flow component [13], 

(tt + ipq 2 x )((l + kl)tt 2 -a 2 ) = K [q x (l + kl)n-2k x a}\ ( f> DW \ 2 , (6) 

K = 2pk 2 y q x k ± _, (7) 

5 



k± _ = kl + kl- q 2 x , (8) 
(3 = 1 + r + rS, (9) 



u + \Tt n gik y 



a = (1 + feju; - fcj, + e n gik y (3. (11) 

Here, Q is the zonal flow real frequency and growth rate and q x is the zonal 
flow wavenumber. In this letter, the third order dispersion relation is the 
starting point for the calculating the back-reaction of zonal flows on the ITG 
mode turbulence. 

To calculate the back-reaction, we start by calculating a small pertur- 
bation to the linear state using Eq. 1 and utilizing the dominant nonlinear 
coupling of the sidebands to the zonal flow to back-react on the turbulence, 
This can be written as, 

f iDW = L Q (j) + —^((fi + (f)* ZF - 4>*_4>zf)- (12) 
e n giT 

Here L is the linear part of the solution, 

lu + (u + aiky)k\ - (1 - (1 + T)e n gi)k y 

~ T^k ' [ } 

It is also assumed that in calculating the non-linearities the linear relation- 
ship between the electric potential perturbation (0) and the ion temperature 
perturbation {%) can be used. Here, it is important to note that the positive 
and negative sideband couples through the non-linear terms only. 

Next, the sideband relationships are considered. The equations for the 
positive and negative sidebands can be written, 

7 -ail 4>DW<PzF /-.,>. 

<p + = '/Wx- (1 + A;i)n + a , (M) 

*- - -Wx- (1 ^ o (15, 



6 



The effects of the ion temperature non-linearity ([</>, Tj\) are neglected, if 
this is taken into account the sidebands become coupled. Moreover, if the 
full system is studied including the ion temperature non-linearity the result 
found here is recovered in an expansion in the zonal flow wave number (q x ) 
and keeping the leading order term. The new real frequency and growth rate 
found from the perturbed linear state is renamed (a)). Utilizing the sideband 
Eq.s 14 and 15 into the perturbation to the linear state Eq. 12 the modified 
temperature perturbation can be found, 

Tiow — Lq4> dw 

~ Kl{ (l + kl)Q + «-(! + m& - a ^ ZF ^ DW ' (16) 
= L cj) DW - NL\4> ZF \ 2 (f) DW , (17) 

K x = (5 2 q x k 2 y ^- (18) 
Eq. 17 is now used together with the linear part of Eq. 2, 



fa ~ l) k y + l u I 
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Tiow — —I — - — r~ 4>dw- (19) 



Combining Eq.s 17 and 19 a modified second order algebraic equation for the 
linear state modified by the zonal flow can be found, 

1 0r 

(1 + k 2 ± )u; 2 + (-(1 - (1 + -^-)e n9i - 

-k 2 ± (ai + ^Te n gi)) + ^\(pzF\ 2 )kyU> + 

7 5 5 5iVL 

{m-- + -{l + r)t n g i + -a l kl + ^\^ ZF \ 2 )k 2 y Tt n g i = 0. (20) 

This equation describes a linear drift wave state modified by the effects of 
zonal flows. Moreover, it is seen that the zonal flow saturation level can 
directly modify the ITG mode stability threshold that may correspond to 
the so-called Dimits shift, observed in numerous non-linear transport simu- 
lations 1 161. 



Ill Saturation levels 



To quantify the effects of the zonal flow suppression an estimate of the zonal 
flow saturation level is needed. The zonal flow saturation may be subject to 
several nonlinear mechanisms in addition to the neoclassic collisional damp- 
ing. In a collisionless plasma, possible nonlinear processes to limit the zonal 
flow instability include the spectral modulation of the ambient turbulence 
due to the flow, tertiary instability like the familiar Kelvin-Helmholtz(KH) 
mode, the nonlinear wave-packet scattering and/or trapping [6], [11]. Here 
to sidestep the complex, specific nonlinear process, an approximate method 
similar to the mixing length theory for estimating the drift wave turbulence 
saturation level, which involves the spectral modulation process of the turbu- 
lence, is applied to roughly evaluate the zonal flow level and the parametric 
dependence on the turbulent pump fluctuation. In previous research the drift 
wave saturation level has been estimated through mixing length theory, 
1 
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This estimate gives a relation between the drift wave electric potential per- 
turbation and the temperature gradient length scale. It is assumed that the 
main source of free energy is the temperature gradient. The zonal flow sat- 
uration level can be estimated through a balance of the vorticity equation. 
The mode coupling saturation level for drift waves is found by balancing the 
linear drift wave growth rate with the main (E x B convection) non-linearity. 
In analogy, the zonal flow saturation level is estimated by balancing the linear 

growth rate with the Reynolds stress term, 
dV\~ 2 dj DW 84> ± 

~W* ZF ~ V ^^^> (22) 

It is argued that the saturation of the zonal flow instability occurs when 
the sidebands become comparable to the pump wave i.e., \4>±\ ~ |0dw|. 
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Hence, 
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Here A is a proportionality constant. From Eq. 23 it is possible to calculate 
the zonal flow saturation level. In the last equality Eq. 21 has been used. 
This saturation level gives some interesting insights, the zonal flow saturation 
level is increased by a streamer like source drift wave (large k y and small k x ) 
similar results were found in Ref. [12] and it is similar to previous results 
found in gyrokinetic theory [27]- [28] . 

To corroborate the analytical estimation of the zonal flow saturation 
level numerical calculations based on the Hasegawa-Mima turbulence system, 
which is extensively employed to investigate the zonal flow generation [6]- [9] 
and also saturation level [6], [9], have been performed for the equation [8]- [9] 



Here 5 = 1(0) for the zonal flow (finite k y ) components, respectively. In the 
calculation, pump energy decreases as the zonal flow instability is excited 
while the total energy is conserved. We change the total pump energy under 
fixing k x (= 0.125 or 0.25) and k y = 1. In the Hasegawa-Mima system a non- 
linear saturation of zonal flows is present and a good agreement with the 
proportionality 4>zf1zf ~ |0rw| 2 in (Eq- 23) is found, as shown in Figure 1. 
Hence, although Eq.(23) is not rigorously justified analytically, the numerical 
calculations show that the relation between the zonal flow saturation level 
and the pump ITG fluctuation is qualitatively correct and is quantitatively 
represented by the proportionality constant, which may involve the spectral 
relation and other physical processes. 




(24) 



using an assumed pump wave of the form 



(ftp = A p sin(k x x) cos(k y y). 



(25) 
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Figure 1: The zonal flow saturation level as a function of drift wave pump 
energy. k x = 0.125 (open circles) or k x = 0.25 (closed circles), k y = 1.0, 
q x = 0.75. 

IV Result and discussion 

The algebraic dispersion relation (Eq. 20), describing the effects of zonal flow 
back-reaction on the background turbulence, is solved in two cases using the 
derived saturation levels for drift waves and zonal flows. 

In Figure 2 the drift wave growth rate as a function of normalized temper- 
ature gradient (rji) is shown with zonal flow wave number (q x ) as a parameter 
for r = 1, k x = k y = 0.3, e n = 1.0, /i = 0. The results are displayed for 
q x = 0.0 (squares) and q x = 0.8 (plus). Figure 2 shows that the zonal flow 
has a significant suppression effect on the drift wave growth rates close to the 
threshold, the growth rates are reduced by > 50% in the range 1.8 < rji < 2.3. 
The critical gradient is also modestly increased by the zonal flow suppression. 

Figure 3 illustrates the effect of zonal flow back-reaction on the critical 
gradient (rj it h) in the e n - rji plane for short wave length zonal flows (q x = 0.8). 
The parameters are r = 1, k x = k y = 0.3, fj, — 0, q x — 0.0 (squares) and 
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Figure 2: The drift wave growth rate and real frequency as a function of 
normalized temperature gradient (77$) is shown with zonal flow wave number 
(q x ) as a parameter for t = 1, k x = k y = 0.3, e n = 1.0, /i = 0. The results 
are displayed for g x = 0.0 (squares) and q x = 0.8 (plus). 

q x = 0.8 (asterisk). It is found that the critical gradient is significantly 
increased for peaked density profiles (small e n ) relevant for edge plasmas. 
For both cases the FLR effects tend to decrease the stability threshold for 
peaked density profiles whereas the non-linear contribution tends to increase 
the threshold. 

At present, the temperature non-linearity has been neglected for simplic- 
ity, if it is incorporated more involved sideband expressions will be found 
where the sidebands couple, however, the effect will be a factor q\ smaller 
than the terms included here. 

The effect of zonal flow suppression is generally manifested in a decreas- 
ing linear growth rate and/or the up-shift of the critical gradient (stability 
threshold). There are other possible stabilization mechanisms and it is inter- 
esting to note that a change in the spectrum due to the back-reaction of the 
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Figure 3: The effect of zonal flow back-reaction on the critical gradient (rjuh) 
in the e n - rji plane for short wave length zonal flows (q x = 0.8) is shown. The 
parameters are t = 1, k x = k y = 0.3, // = 0, q x — 0.0 (squares) and q x = 0.8 
(asterisk). 

zonal flow may also change the diffusion, moreover, it has been shown that 
the spectrum as well as the drift wave growth rate are quite prone to plasma 
shaping effects [13]- [H] which is out of the scope of the present study. 

V Conclusion 

In summary, the present study aims to clarify the explicit effects of zonal 
flow suppression on drift waves based on a quasilinear theory. The short 
wave length zonal flow suppression of drift wave turbulence is significant in 
particular close to the stability threshold (rjuh)- In the regime above the linear 
threshold the drift wave growth rate is significantly reduced indicating a 
significantly reduced transport. To explain the total quenching of turbulence 
in the wider Dimits shift region a more complete model is needed including 
more physics as well as plasma shaping effects. Numerical calculations of 
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realistic ITG turbulence with focus on the zonal flow saturation level is also 
planned. These will be treated in future papers. 
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